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ABSTRACT 

We analyzed the radial surface brightness profile of the spiral galaxy NGC 7793 using HST/ ACS images from 
the GHOSTS survey and a new HSTfW¥C3 image across the disk break. We used the photometry of resolved 
stars to select distinct populations covering a wide range of stellar ages. We found breaks in the radial profiles 
of all stellar populations at 280''(~5.1 kpc). Beyond this disk break, the profiles become steeper for younger 
populations. This same trend is seen in numerical simulations where the outer disk is formed almost entirely 
by radial migration. We also found that the older stars of NGC 7793 extend significantly farther than the 
underlying Hi disk. They are thus unlikely to have formed entirely at their current radii, unless the gas disk was 
substantially larger in the past. These observations thus provide evidence for substantial stellar radial migration 
in late-type disks. 

Keywords: galaxies: evolution — galaxies: spiral — galaxies: stellar content — galaxies: structure — tech- 
niques: photometric 



1. introduction 

Studies of the chemical composition of galaxy disks typ- 
ically assume that the orbital radii of stars are essentially 
static. Only a small diffusion in radial distance, increasing 
with age, is expected due to dynamical heating (e.g., |WieTen 
p^977( pehnen & Binneynl998| . In this scenario, the metal- 
licity of stars of similar ige and galactic distance should be 
comparable, as they are bom from the remnants of previous 
stellar generations that orbited at the same radius. However, 
observations of local F and G dwarfs have found little evi- 
dence for such coexistent evolution, instea d revealing a large 
scatter in the age-metallicity distribution ( [Edvardsson et al. 
1993| ). Theoretical work has suggested that this scatter may 



be due to resonant interactions between stars and spiral arms, 
which enable stars to migrate across many kpc ( Sell wood & 
Binney'2002'). Indeed, using a high-resolution hydrodynam- 



ical simulation, |Roskar et al.| (E008a ) found that such migra- 
tion naturally led to the high degree of variation observed in 
the age-metallicity relation of local Milky Way (MW) stars. 
These events, in which the guiding center of the star shifts, 
differ from general diffusion where the epicyclic motion in- 
creases by only a few kpc due to an increase in velocity dis- 
persion. 

Testing such models of stellar migration may be possible 
using studies of stellar populations beyond the classical disk 
break of late-type galaxies. The cause of this break is still de- 
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bated (see the review by [Trujillo & Pohlen||2QQ5| ), however 
the change in surface brightness is often ass ociated with a 
threshold in the star formation rate (SFR; e.g., "van der Kruit 
& Searle 198T1 [Kenhicutt 1989; Schaye 2004; Elmegreen & 
Hunter 2006 ). Using a study of 85 late-type galaxies, Pohlen 
& Trujillo (2006 ) found that the surface brightness profile o 



approximately 60% of disk galaxies steepens beyond the disk 
break. Hence, the outskirts of disk galaxies are often associ- 
ated with a dearth of star form ation. Using A^-body simula- 
tions of disk galaxy formation, |Roskar et al.| ([2008b ) found 
that if star formation is sufficiently suppressed beyond the 
break, then that region may form almost entirely from stars 
migrated from the inner disk. This would leave a clear sig- 
nal on radial surface brightness profiles, with older stars (i.e., 
those that have been subject to scattering for longer) exhibit- 
ing shallower profiles. 

Observati ons from the GHOSTSTj survey ( Radburn-Smith 
et al.||20lT[ hereafter R-Sll), which have targeted the out- 



skirts of several nearby disk galaxies using the Hubble Space 
Telescope Advanced Camera for Surveys (HST /ACS), thus 
offer a significant test of radial migration. GHOSTS obser- 
vations of NGC 4244 revealed the location of the disk break 
to be independent of stellar age, thus constr aining several 
mechanisms for the formation of this feature ( |de Jong et al. 
2007| ). However, lack of data beyond the break prevented fur- 
ther analysis. In this paper, we focus on the spiral galaxy 
NGC 7793. This flocculent galaxy is (1) relatively nearby 
(3.7±0.1Mpc, R-Sll), (2) situated at high Galactic latitude 
(-77.2°) and hence relatively uncontaminated by foreground 
stars, and (3) an isolated loosely bound member of the Sculp- 
tor Gr oup. NGC 7793 was recently studied by |Vlajic et"aL 
( |2011| ), who found no sign of a break in the red giant branch 
(RGB) population, but did find evidence for an upturn in 
the radial metallicity distribution beyond the disk break. In 
this paper, we expand on this work by using data from the 
GHOSTS survey together with a new HST Wide Field Cam- 
era 3 (WFC3) image of the inner disk to measure fainter stars 
and to recover stellar densities in this crowded region. We 
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Figure 1. Location of the observations and resulting color magnitude diagrams (CMDs). Radial surface brightness profiles are taken from the shaded wedge 
(with a de-projected opening angle of 25°) in panel (a). The thick red line marks the disk break. The blue dashed region denotes the area covered by the 
Spitzer/IRAC 8 fim observation in the top right inset, which highlights the spiral structure that may be responsible for stellar migration. Panels (b)-(d) plot 
the individual CMDs, with the regions used to extract the different stellar populations, as described in the text, highlighted in panel (b). Gray and black points 
indicate different cuts on the signal-to-noise ratio of 3.5 and 5.0, respectively, with contours drawn at densities of 50, 100, 150, and 250 stars per 0.01 mag^. The 
black arrows indicate the direction of reddening due to foreground extinction, which has been removed from the data. The shaded ellipses indicate photometric 
errors as reported by dolphot, and the blue dashed lines mark the 50% completeness limits calculated from artificial star tests. 



compare these observations with predictions from the radial 
migration simulations of Roskar et al. (2008b) and discuss al- 
ternate explanations. 

2. //^r ACSAVFC3 PHOTOMETRY 

Figure [T| plots the location of the two GHOSTS fields dis- 
cussed in this paper (Fields 2 and 3). Each field was observed 
with dithered 740 s exposures in both the F606W and F814W 
filters. R-Sll described the data reduction and photometry 
of the resolved stars in these fields using the DOLPHOT code 
( |Dolphinl[2QQQ| ). A further field (Field 1) was observed with 
the WFC3 camera to characterize the disk truncation (HST 
GO program 12196). This observation was made using the 
same filters as the ACS fields, with dithered exposures of 
980 s in F606W and 1488 s in F814W. We use the same 
DOLPHOT parameters and the same crowded-field stellar se- 
lection criteria for the WFC3 field as used by R-Sll for 
the ACS photometry. Despite minor differences between the 
cameras, namely, pixel scale and quantum efficiency in the 
F814W filter, we found little discrepancy in the resulting pho- 
tometry between the cameras in overlapping regions. 

The resulting color magnitude diagrams (CMDs) for each 
of the fields are presented in Figure [T] By analyzing syn- 
thetic CMDs generated with a constant SFR of 0.01 



yr~\ an appropriate metallicity of [Fe/H]= -1, and photomet- 
ric errors matching the observations, we identified three re- 
gions of the CMD corresponding to three discrete age ranges. 
The area marked as "Young" consists of bright main se- 
quence and helium-burning (HeB) stars, which are collec- 
tively 10-100 Myr old; the "Intermediate" region is composed 
of 100 Myr -1 Gyr old HeB stars; and the "Old" region con- 
sists of 1 - 10 Gyr old RGB stars. We subsequently use these 
regions, which suffer little cross contamination, to select stel- 
lar populations by age. 

R-Sll generated artificial stars for the ACS fields, which 
were subsequently used to assess errors and correct for com- 
pleteness due to stellar crowding. We generated an additional 
1.8 million artificial stars for the WFC3 image. Figure [T] plots 
the 50% recovery limit of these artificial stars, which we set 
as our limiting magnitude. 

We found good agreement between detection rates and pho- 
tometry in the overlapping region of Fields 1 and 2. The 
WFC3 observation from Field 1 yielded 3474 detections 
brighter than F814W = 25.75 mag after culling, while ACS 
(Field 2) found 3631, despite the shorter exposure time. Ap- 
proximately 82% of the ACS detections were coincident be- 
tween fields. After correcting for completeness with the ar- 
tificial star tests, this increased to 95%. The remainder can 
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Figure 2. Radial surface brightness profiles of "Old" (red), "Intermediate" aged (yellow), and "Young" (blue) stellar populations. The left panel plots number 
counts in NGC 7793 taken from the shaded region of Figurel^ panel (a). The middle and right panels, respectively, plot surface densities extracted from similar 
regions of the Roskar et al. ( 2008b I simulation of a Milky Way (MW) sized galaxy and the new simulation of NGC 7793 described in the text. The NGC 7793 
observations have been linearly offset such that the "Old" population matches a Spitzer 3.6 /im surface brightness profile, indicated by a solid red line. For clarity, 
the "Intermediate" and "Young" population are arbitrarily offset by 5 and 12.5 mag, respectively, in the observations, and by 2 and 4 dex in the simulations. 
All profiles are scaled to the approximate location of the "Young" break (280^^in NGC 7793, 10.2 kpc in the MW simulation, and 7.7 kpc in the NGC 7793 
simulation). Arrows indicate where the Poisson uncertainty may move the data point below the background level. Dashed lines indicate exponential fits to the 
inner and outer profiles, with the corresponding normalized exponential-disk scale lengths given by Ri and Ro, respectively. The shaded regions in the left panel 
indicate the range of possible profiles due to different treatments of background contamination as described in the text. Faint solid lines in the simulation panels 
plot the radial surface brightness profiles using the formation radius of the stars, i.e., removing the effects of stellar migration. 



be attributed to incompleteness at the magnitude limit of the 
survey. The median photometric differences between coinci- 
dent detections are -0.05 and 0.03 mag in F606W and F814W, 
respectively. The mean error of these detections, when com- 
bined in quadrature, is ±0.10 in both filters. Hence, although 
the offsets are systematic, they are significantly less than the 
DOLPHOT uncertainties. We thus found our photometry from 
the two different cameras to be consistent and complete to 
F814W = 25.75 mag. 

3. RADIAL SURFACE BRIGHTNESS PROFILES 

In the left panel of Figure |2] we plot the radial star counts 
of NGC 7793 within the shaoed region of panel (a) of Fig- 
ure [T] for our key stellar populations. These profiles are all 
normalized to the radial break observed in the youngest pop- 
ulation. In addition to masking foreground stars and back- 
ground galaxies, a region of significant star formation, located 
at 23^58^13^1, -32°36'37;'3 was masked out of Fields 1 and 
2 to provide a better estimate of the mean surface brightness. 
We assigned an equivalent surface brightness to the star count 
profiles by multiplying the logarithm of the older RGB counts 
by 2.5 and linearly offsetting the result to match a Spitzer In- 
frared Array Camera (IRAC) 3.6 fim surface brightness pro- 
file extracted from the same region. Consequently, we find 
the star count profiles reach a limiting magnitude equivalent 
to ~ 32 AB mag arcsec"^ at 3.6 /im. 

Each of the stellar population bins denoted in panel (b) of 
Figure [T] are substantially brighter than the completeness lim- 
its. They are thus less susceptible to spurious detections from 
cosmic rays and CCD noise that become significant at a lower 
signal-to-noise ratio. Instead, contamination is expected from 
unresolved background galaxies, which far outnumber any 
foreground MW stars. To account for this, we followed the 
same analysis used by R-Sl 1 and ran the photometry code on 
ACS images of high-z sources with exposure lengths longer 



than our observations. As these fields are free of any resolv- 
able stars outside of the MW, any detections will thus be con- 
taminants, primarily from background galaxies. For each of 
the denoted stellar populations in panel (b) of Figure [T] we 
measured contamination from the number density of remain- 
ing detections in these "empty" archival fields after imposing 
the selection criteria used by R-Sl 1. For the "Old", "Interme- 
diate", and "Young" bins we recovered 4.4 x 10""^, 4.0 x 10""^, 
and 3.3 x 10""^ detections per arcsec^, respectively. Each of 
these contaminant measurements are lower than any of the 
detection rates in the corresponding NGC 7793 populations. 
Averaging the outermost radial bins of each population in 
NGC 7793, we, respectively, measured 4.0 x 10"^ 4.5 x 10"^, 
and 1.6 X 10 ~^ detections per arcsec^. As discussed in lElHs &J 



Bland-Hawt horn (2007), these can be treated as upper limits 
on the true level of contaminants in each bin. Moreover, us- 



ing the Galaxy Count software described in |Ellis & Bland- 
|Hawthorn| ( |2007 |, we can infer uncertainties on these values 
due to cosmic variance and Poisson noise. These equate to 
upper limits on the total number of background galaxies in an 
ACSAVFC field of 163 ± 14, 18 ± 4, and 65 ± 8 for the re- 
spective age bins. These uncertainties are incorporated into 
the final error we assign to the profiles. The lower bound on 
the profiles due to removing this higher level of contamina- 
tion, as well as an upper bound resulting from not subtracting 
any contaminants, are indicated by the shaded regions in the 
left panel of Figure [2] 

We are able to probe the outskirts of NGC 7793 much 
deeper than previous integrated light studies, whose profiles 
do not cov er these radii o r reach this surface brightness limit 
(e.g., Munoz -Mateos et a l. 2009; Vlajic et al. 2011 ). We con- 
sequently foundThat the disk break, ^break, occurs at approxi- 
mately the same location regardless of the age of the underly- 
ing stars, in agreement with observations of the edge-on disk 
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Table 1 

Exponential Outer-disk Fits 



D^ta Set ^old /-^break ^int / -^break %oung / ^break 

NGC 7793 Observations 

Raw 0.147±0.007 0.107±0.016 0.078±0.014 

Corrected 0.145±0.007 0.098±0.016 0.075±0.014 

Overcorrected 0.134±0.016 0.080±0.017 0.063±0.022 

Simulations 



Roskar et al. 



1 2008b 



INUC / /93 Model 



0.134±0.008 
0.125±0.012 



0.072±0.007 
0.101±0.007 



0.053±0.006 
0.094±0.007 



Roskar et al, 



Simulations without Migration 



2008b 



INUC / /93 Model 



0.025±0.018 
0.054±0.019 



0.018±0.011 
0.032±0.011 



0.024±0.014 
0.025±0.013 



Note. — The "Raw" NGC 7793 observations refers to profiles uncor- 
rected for contamination from background galaxies. For the "Overcor- 
rected" profiles, an upper limit on the number of the contaminants is in- 
ferred from the average density of detections in the outermost radial bins 
of each stellar population in the raw data. The exponential-disk fits to the 
"Simulations without Migration" use the formation radii of stars born in the 
outer-disk region. 



galaxy NGC 4244 ( |de Jong et aH2Q07| ). This bre ak occurs at 
28C(5.1 kpc), internal to the regions covered by Vlajic et al. 
([2011) and external t o the integr ated light profiles of Munoz-^ 



Mateos et aLlpQQ9| ). |de Jong et a l. ( 2007 ) argued that a con 



stant break location calls for a dynamical interpretation, given 
that a break radius due solely to a star formation threshold is 
unlikely to remain at the same radius over the lifetime of the 
galaxy, due to disk growth. We note that the break is weakest 
in the oldest population, which also shows a possible second 
break to a shallower profile around ~ 1.8 x T^break- This may 
correspond to a transition in the number counts from a disk- 
dominated population to one composed primarily of stellar 
halo stars (see Section [STT] ). 

Although the disk break appears at a constant radius for 
all stellar ages, the profiles in Figured] suggest that the slope 
of the profile does depend on age, botn within and beyond the 
break. We characterize this slope by fitting the inner and outer 
profiles to an exponential surface brightness function of the 
form /i(r) = j^oe'^^^^ where /io is the central surface bright- 
ness and Rd is the exponential-disk scale length. The resulting 
fits, shown in Figure [2] and Table [T] reveal a smooth increase 
in the steepness of the outer profile with decreasing age, to- 
gether with a corresponding flattening of the inner profile. 

To assess the uncertainty in the profiles we use Monte Carlo 
realizations of the data. As summarized in Table [T] the outer- 
disk parameters of the "Old" population are robust (~5% un- 
certainty), but low number counts in the "Intermediate" and 
"Young" outer-disk profiles leads to uncertainties of ~16% 
and ~19%, respectively. Accounting for different levels of 
contamination has little effect on the fits. 

4. COMPARISON WITH SIMULATIONS 



In isolated simulations of a MW-type galaxy, |Roskar et~aL 
( [2008 b ) found significant migration of stars over several kpc, 
such that ~90% of the outer disk was predominantly formed 
from migrated stars. Older populations exhibited larger scale 
lengths than younger populations, because the former were 
subject to migration for l onger. We can compa re our mea- 
sured profiles to those of [Roskar e^al.| ( |2008b| ), which are 



plotted in the middle panel of Figure [2|f or the same age ranges 
used in our CMD analysis. These profiles were extracted from 
the simulation at the same inclination as NGC 7793 (~53°), 



and for comparison with the observations are similarly scaled 
to the break radius of the youngest stellar population. The 
behavior of the profiles in the simulation is qualitatively sim- 
ilar to what is observed in the data. However, the simula- 
tion is not necessarily the best analog for NGC 7793. The 
mass of the simulation contained within 10 kpc is approx- 
imatel y 1.2 x IO^^Mq, larger th an the mass of NGC 7793, 



which Carignan & Puche| ( |1990| ) measure within 7.35 kpc as 
1.5 X 1O^"M0. Similarly, the radius of the disk break is also 
larger in the simulation (10.2 kpc cf. 5. 1 kpc in NGC 7793). 

4.1. Modeling NGC 7793 
A potentially better analog to NGC 7793 i s the equivalent 



simulation generated to model NGC 300 in Gogarten et al. 
( |2010| ). This galaxy is another member of the Sculptor Group 
with mass similar to NGC 7793, although the star formation 
history (SFH) of the system and the lack of a disk break in- 
dicate a quiescent formation. Gogarten et al. (2010) found 
little evidence of migration in HST observations of NGC 300, 



and le ss migration in the NGC 300 simulation than the |Roskar 
et al.|p008b ) model. This led the authors to suggest that both 



mass and environment may be factors in radial migration. As 
the formation of transient spiral structur e is linked to gas ac- 
cretion and star formation ( [Sellwood & Carlberg 19 84), such 
a dependence on enviro nment is to be expec ted. This was 
further demonstrated by Gogarten et al. ( |2010| ), who note that 
the comparably massive M 33, whic h is interacting with M 3 1 , 
does show signs of radial diffusion ('Willia ms et al.|2009j ). 

We extend this analysis by observing that the ^-band abso- 
lute surface brightness of NGC 7793 is on average ^l.Oma^ 
arcsec"^ brighter than NGC 300 (see panel (a) of Figure |3]). 
The surface density of the NGC 300 disk may thus be too 
low to support the transient spiral structures required to drive 
substantial migration. To better match the properties of 
NGC 7793 we resimulated the NGC 300 model with a smaller 
spin parameter to yield a disk with a smaller scale length 
(Figure [3] panel (b)), larger rotational velocity (Figure [3j 
panel (d)), and higher surface density within the central 4 kpc. 
These changes resulted in a more concentrated system, with 
1.7 X 10^^ M0 contained within the disk break (located at 
7.7 kpc). The stellar disk in this model was also subject to sig- 
nificant radial migration, leading to a s imilar increase in disk 



scale length with stellar age as in the [Roskar et al.| ( [2008b 
model (see right panel of Figure p| and Table [jj). The effect is 
however, weaker. We partly attrioute this difference to a rela- 
tive increase in the level of in situ formation beyond the disk 
break in the NGC 7793 model. This is manifest in Figure [2[ 
and Table [T] by exponential-disk fits to the formation radii of 
stars located beyond the break. Specifically, the normalized 
scale lengths are greater in the NGC 7793 model than in the 
MW simulation. The NGC 300 model generated by Gog arten[ 
[et al. ( 2010 ) shows a further increase in this in situ population, 
which resulted in an even softer disk break. 

Figure [4] compares the cumulative SFH of the main disk 
of NGC 7793 to the ,Roskar et al.| ( |2008b) MW simulation 
and the new NGC 7793 simulation, both of which show signs 
of radial migration. To extract the SFH from the HST / ACS 
photometry, we used the software package match (Dol phin[ 
2002[ ). This program uses a maximum likelihood analysis 
to tit the observed CMD with stars generated from the stel- 
lar evolution models of IGirardi et al.| (120021 with updated 



lar evolution models of Girardi et al. ( 2002 with updated 
AGB models from Marigo et al." 200 8 "and Girar di et al.|2QTQl 
and assuming a ,Salpeter .1955. IMF and a binary fraction of 
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Figure 3. Radial surface brightness profiles (top panels) and Hi rotational velocity profiles (bottom panels) from observations of NGC 7793 and NGC 300 (left 
panels) and from simulations of NGC 7793 and NGC 300 (right panels). The NGC 7793 simulation is described in the text while the NGC 300 simulation is 
taken from Gogarten et al. (2010). The K-bmd surface brightnesses plotted in panel (a) were measured by Munoz-Mateos et al. ( 2007) from the Two Micron 
All Sky Survey Extended Source Catalog ( Jarrett et al. 2000). The shaded regions indicate the la uncertainty in the photometry from variations in intensity 
within each isophote and errors in sky subtraction. The Hi rotation curves were measured by Carignan & Puche (1990) and Puche et al. (1990) from Very Large 
Array observations of NGC 7793 and NGC 300, respectively. The shaded regions indicate the velocity uncertainties derived from the largest difference between 
the average radial rotation and the rotation curve measured on each side of the galaxy. Data are shown to the disk break in NGC 7793 for both the observations 
(~ 5.1 kpc) and simulations (~ 7.7 kpc). 



40%). The artificial CMD was calculated using a 36x25 grid 
in log(age) - [Z] space with bin widths of 0.1 dex. We left the 
distance modulus and extinction free to be fitted by MATCH, 
which converged on a distance modulus of m-M = 27.83 and 
an extinction of Ay = 0.05. These are consistent with the lit- 
erature values of m-M = 27.87 ±0.08 (R-Sl 1) and Ay = 0.06 
( [Schlegel et al.||1998] ). We constrained the chemical enrich- 
ment history to an increasing metallicity. 

The overall uncertainty in the SFH, denoted in Figure [4] by 
the shaded region, is a combination of both random and sys- 
tematic uncertainties. The former are measured from the la 
spread in 300 Monte Carlo realizations of the best-fit solution. 
The systematic uncertainties in the isochrones are calculated 
by varying the Mboi and T^ff values of the underlying models 
(for details, see Dolphinj2012 ). 

As the simulations are only evolved for 10 Gyr, we tie their 
initial cumulative star formation to the equivalent value mea- 
sured in NGC 7793 10 Gyr ago. The resulting SFHs of the 
models show a relatively smooth build up of mass. This is to 
be expected, as the simulations are not set in a full cosmolog- 
ical environment. The SFH of NGC 7793, as reconstructed 
from the observations, is consistent with a similar steady SFR 
for the life of the system. 

Even though the SFHs of the two simulations are similar, 
they exhibit markedly different inner-disk profiles (see Fig- 
ure [2]). This suggests dynamical processes dominate the con- 
struction of the disk profil es, wheth er internally by secular 
processes (as in Roskar et al.|2008b| ) or by externally driven 
merger processes related to the formation of the host halo 



(e.g., the two controlled simulations presented by |Purcell et al. 
201 1| ). Outer-disk profiles are likely similarly affected, and in 
the case of the simulations presented here, they predominantly 
arise from radial diffusion. Figure [5]presents a direct compar- 
ison of these normalized outer-disk gradients as a function of 
age. There is remarkable agreement between the simulations 
and the NGC 7793 observations, which thus strongly attests 
to significant migration in NGC 7793. 

4.2. Alternate Simulations 



'Sanchez- Blazquez et al.| ( |2009| ) have also addressed radial 
migration in their cosmological simulations. They argue that 
the upturn in the age profile beyond the disk break is prin- 
cipally due to an extended and constant SFR in the outer 
disk. However, they do find evidence for some migration, and 
moreover have studied the outer regions of their cosmological 
simulations without the effects of such migration (see their 
Figure 17). Specifically, they found a break still formed in the 
radial surface brightness profiles, but without migration, the 
decrease in gradient of the outer-disk profile with increasing 
stellar age was lost. 

Stellar radial migration may also be driven by other mech- 
anisms, such as resonances between spiral structures and cen- 
tral bar^ (Minchev & Famaey 2010 ), or interactions with 
satellites (Quillen et al. 2009; Bird et al. 2012). These pro- 
cesses may operate in addition to the resonant scattering from 

^ ^ Although we note that resonances with the central bar are likely to have 
less influence on the outer disk (Debattista et al. 2006 Minchev et al. 20T2| . 
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Figure 4. Comparison of the cumulative star formation history of NGC 7793 
with the MW and NGC 7793 simulations as a function of age. The SFH 
for NGC 7793 was generated using MATCH, with the shaded region denoting 
the uncertainty in the measurement. As the simulations are evolved for only 
10 Gyr, we fix their initial cumulative star formation to that of NGC 7793 
10 Gyr ago. A zoom-in of the last 600 Myr is shown as an inset. The isolated 
simulations show a smooth buildup of stellar mass over the lifetime of the 
systems, which is consistent with the SFH reconstructed from observations 
of NGC 7793. 



transient spiral arms, and so serve to amplify the effect. How- 
ever, as NGC 7793 is not host to a si gnificant bar, and as such 
features are believed to be robust (jShen & Sellwood||2Q04t 
[Athanassoula et al.|[2005l [Debattista et al. || 2Q06| ), the disk is 
unlikely to have been affe cted by a central bar. Conversely, 
as detailed in Section [53) NGC 7793 may have undergone a 
recent merger. Such an event could enhance spiral structure 
in the disk and so increase the level of stellar migration. 



5. ALTERNATE EXPLANATIONS 

The increase in the outer-disk scale length with stellar age, 
as seen in Figure [5] corresponds to an increase in the aver- 
age age of the outer disk with radius. This trend, although 
not ubiquitous, has been confirmed using spectroscopic stel - 
lar ages in several nearby disk galaxies ( Yoachim et al.|2Q12] ). 
To assess whether radial migration is the only viable mecha- 
nism for setting this trend, we consider alternate explanations 
below. 

5.1. Thick Disk and Halo Contamination 

The thick disk typically extends farther radially than the 
thin disk by a factor of 1.25, and in low-mass galaxi es can ac- 
count for n early hal f the total ma ss and luminosity ( [Yoachim 
|& Dalcan ton 2006 ; [ Comeron et al. 2011). As the thick and 
thin disks can only be separated in high inclination systems, 
a significant portion of the older outer-disk stars may actually 
lie in an extended thick disk component. The smooth change 
in gradient of the outer radial profile with age argues against 
a simple two component model. However, if true, the coinci- 
dent break radius of both old and young populations suggests 
that the mechanism propagating the break, and thus setting 
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Figure 5. Comparison of the normalized exponential-disk scale lengths mea- 
sured in Figure[2]for the outer disk of NGC 7793, the Roskar et al. ( 20()8bJ 
MW simulation, and the new NGC 7793 simulation. The width of each box 
represents the range of ages covered by each population, while the height en- 
codes the Icr uncertainty in the exponential-disk fit as calculated from Monte 
Carlo resamples. The width of the NGC 7793 simulation boxes are slightly 
reduced for clarity. Both observations and simulations argue for a smooth 
increase in the scale length of the outer-disk profile with stellar age. 

the outer-disk gradients, applies equally to both components. 
Furthermore, the thick dis k itself m ay have form ed from a mi- 
grated population ( Schonrich & Binney 2009 ; Loebman et al. 
[201 1| . Hence, contamination from a thick disk is unlikely to 
affect the inference of migration. 

In the "Old" RGB population of NGC 7793 plotted in Fig- 
ure [2j a gradual departure from the outer exponential is seen 
below 29 mag arcsec"^. This may indicate a transition to 
either a thick disk or old stellar halo component. At r = 
1 1 kpc (2.2 X 7?break), this extra component accounts for some 
4.5 X 10~^ stars arcsec"^. This is approximately a factor of 10 
greater than the number of predicted contaminants, and simi- 
lar to the surface densities attributed to the stellar halos of iso- 



lated si mulations of disk-satellite interactions from Pur cell,| 
[Bullock & Kazantzidis (2010). However, such low numbers 



will have minimal impact on the exponential fits summarized 
in Figure [5] 

5.2. In Situ Formation 

Star formation in o uter disks has been characterized by 
Bigiel et al. (2010a|b), who find evidence for spatial corre- 



lation between far-ultraviolet (FUV) emission, which traces 
recent star formation, and emission from atomic hydrogen gas 
(Hi), which compris es most of th e interstellar medium in the 
outer disk. Although [Bigiel et aL] ( f2010b) find the overall star 
formation efficiency to be extremely low m these regions, they 
observe that this efficiency scales with Hi gas density. Hence, 
as the gas disk of NGC 7793 extends beyond the disk break, 
we expect some of the stars in this region to have formed in 
situ. This may partly account for the vertical offset seen in 
Figur e [5] between the NGC 7793 observations and the [Roskar 
let al.| ([2Q08b) simulation. 

To explore this possibility. Figure [6| plots the surface den- 
sities of Hi gas mass from the THINGS survey ( Walter et al. 
[2008] ), and the SFR inferred from the GALEX Ultraviolet At- 
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Figure 6. Top: the Hi gas surface density in NGC 7793 from the THINGS 
survey. Middle: the star counts corresponding to the old RGB (red) and 
young MS and HeB (blue) profiles shown in Figurej2] Bottom: the star 
formation rate density from GALEX FUV emission. Tne old RGB stars are 
found to extend much farther than the Hi disk. 



las of Nearby Galaxies (Gi l de Paz eFaLl[2007l ). The FUV 
emission does indeed extend slightly beyond the truncation. 
However, using MATCH with this inferred SFR to generate 
artificial CMDs underestimates the number of "Young" stars 
observed by almost an order of magnitude (although this cal- 
culation does carry significant uncertainties). The stars may 
instead have formed in previous episodes of formation not 
probed by the FUV flux. In deed, the gas disk e xtends further 
stiU, although, as noted by |Vlajic et al.|pon ), the Hi disk 
of NGC 7793 is unusually small. As shown m Figure [6] the 
radial extent of the gas is significantly smaller than the RGB 
distribution. Thus, the bulk of the old RGB stars were not 
likely formed in situ from the Hi gas reservoir now present in 
the region, unless the Hi disk was substantially larger in the 
past. 

5.3. Accretion/Harassment 

Many interacting systems are evident in the Local Group. 
Our own galaxy is tidally disrupting the Small and Large 
Magellanic Clouds, produci ng a number of phenomena such 
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IQuillen et aLr2 009 ; Bird et al. 2012[ |Kazantzidis et al. 



2008). In the optical, NGC 7793 does not appear to be un- 
aergoing such a major merger, which would likely disrupt the 
outer-disk profiles. However, the Hi disk of this relatively 
isolated galaxy is warped, and studies have revealed non- 
circular motions on the n orth side, indicative of harassment 
( [Carignan & Puche1|1990| ), or of an advanced stage of merg- 
ing (e.g., |Sancisi|1999| ). Such minor interactions may explain 



the lack of gas in the outskirt s, and may lead to an extended 



and heated componen t (Pu rcell, Bullock & Kazantzidis|2010 



Kazantzidis et al.|2QQ9| ). Moreover, such interactions are also 
likely to induce the transient spiral featu res that drive stellar 
radial mig ration (Pu rcell, Kazantzidis & Bulloc k||2009[ |Pur ' 



|ceU et al.||2011| [Struck, Dobbs & Hwang((20ir 
transient spiral arms can last for a couple of Gyr 



As these 
(and poten- 



as the Magellanic Stream ( Mathewson et aL_ 
[2005^ . Such interactions strongly affect tl 
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tially much longer in star-forming galaxies due to swing am- 
plification processes involving giant molecular clouds, e.g., 
|D' Qnghia et al.|20r2| ), the level of stellar radial migration we 
infer in NGC 7793 may be a direct consequence of a com- 
pleted past merger. 

6. SUMMARY 

We have measured the radial surface brightness profiles of 
three distinct stellar populations across both the inner and 
outer disk of NGC 7793. Breaks in the radial profiles are 
found at the same location regardless of age. However, older 
stars show a steeper profile internal to the break and a shal- 
lower profile beyond the break in comparison to younger 
stars. The observed smooth increase in the outer-disk scale 
length with age suggests a formation mechanism that contin- 
ually operates over a significant fraction of the lifetime of the 
galaxy. 

We have compared the surface brightness profiles to /"/-body 
simulations that exhibit significant radial stellar diffusion, in- 
cluding a new simulation of a system comparable in mass to 
NGC 7793. The gradients of the outer-disk profiles in these 
simulations are set almost entirely by migration due to reso- 
nant scattering. The new simulation shows a slightly higher 
level of in situ formation beyond the disk break relative to 
Roskar et al . ( 2008b ), which corresponds to a flattening of the 
profiles. However, the gradients overall match extremely well 
to the outer disk of NGC 7793. 

We thus argue that the observations presented here are in- 
dicative of high levels of stellar radial migration in NGC 7793. 
This inference carries significant ramifications for studies of 
stellar populations in disk galaxies, including the MW. 

Support for this work was provided by NASA through 
grants GO-10889 and GO-12196 from the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Incorporated, under 
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